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ABSTRACT 

An experimental investigation was performed to determine the effect 
of pressure gradient on the heat transfer to space shuttle reusable 
surface insulation (RSI) tile array gaps under thick, turbulent bound- 
ary layer conditions. Heat transfer and pressure measurements were 
obtained on a curved array of full-scale simulated RSI tiles in a tunnel 
wall boundary layer at a nominal freestream Mach number of 10.3 and 
freestrean unit Reynolds numbers of 1.6, 3.3, and 6.1 x 10^ per meter. 
Transverse pressure gradients were induced over the model surface by 
rotating the curved array with respect to the flow. Definition of the 
tunnel wall boundary layer flow was obtained by measurement of boundary 
layer pitot pressure profiles, and flat plate wall pressure and heat 
transfer. 

Flat plate wall heat trai. .er data were correlated and a method 
was derived for prediction of smooth, curved array heat transfer in the 
highly three-dimensional tunnel wall boundary layer flow. Simulation 
of full-scale space shuttle vehicle pressure gradienx levels was assessed. 
No systematic effect of pressure gradient on RSI tile array gap heat 


transfer was observed- 
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III. INTRODUCTION 


The space shuttle orbiter is being designed for an operational 
life in excess of 100 reentry missions with minimal required refurbish- 
ment between those missions. In order to meet this requirement, the 
vehicle thermal protection system (TPS) will be a surface covering of a 
non-metallic , low-density refractory oxide. This material, referred 
to as Reusable Surface Insulation (RSI), is capable of withstanding 
without degradation repeated e> osure to the harsh reentry environment, 
while insulating the cold structure of the vehicle from surface tempera- 
tures in excess of 1500° K. The material will be attached to the vehicle 
surface in a "brick-like 1 ' array of square tiles (15.24 x 15.24 cm) which 
vary in thickness from approximately 1-10 centimeters according to the 
intensity of the local heating environment. Small gaps between tiles 
will allow for thermal expansion of the tile material. Interference 
heat transfer to the tile exterior and gap wall surfaces is of major 
concern to the TPS designer. The presence of the gaps may result in 
increased boundary-layer turbulence and attendant increased surface 
heating. Heating levels within the gaps, which would be expected to be 
substantially lower than surface values, may be similarly severe due 
to flow rea ! tachment phenomena. In addition, radiation blockage within 
the gaps may produce extreme temperatures even at low heating levels , 
and the shortened heat paths may result in excessive bond-line 


temperatures . 
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Effective design of the TPS requires a sound knowledge of the 
aerodynamic heating environment to which the RSI tiles will be sub- 
jected. This knowledge must include an accurate definition of the heat 
transfer distribution about a tile and a good understanding of how this 
distribution is affected by boundary layer state ( laminar /turbulent ) , 
boundary layer thickness, flow angularity, gap width, tile edge radius, 
tile stacking arrangement, e.nd other parameters. 

Present day understanding of heat transfer in gap or cavity flows 
has resulted from past studies of the general problem of heat transfer 
in regions of separated flow. Chapman (ref. 1) attacked the problem 
of a purely laminar separated flow region theoretically. His analy- 
sis, which assumed a boundary layer of zero thickness * c the separation 
point, predicted that the average heating to the separated region was 
only 56 percent of that for attached flow under similar conditions. 
Chapman’s theory was not capable of predicting heat transfer distribu- 
tions. Larson (ri ■ p . 2) experimentally verified Chapman's result and 
found a similar reduction in average heat transfer of about 60 percent 
for separated turbulent boundary layers . 

Charwat, et al (refs. 3 and 4) made extensive measurements of the 
pressure and heat transfer distributions in cavities under turbulent 
boundary layer conditions at subsonic and low supersonic Mach numbers. 
They were able to identify a critical cavity width-to-depth ratio 
which "separates" the cavity flow into two distinct types. When this 
width-to-depth ratio is exceeded, the flow separates at the forward 
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cavity wall and reattaches to the cavity floor as a rear-facing step, 
then separates again ahead of the upstream-facing cavity val^.. This 
type of flow is termed "closed" cavity flow. When the width-to-depth 
ratio falls below the critical value, the . jundary flow bridges the 
cavity entirely, reattaching at the upstream- facing vail, and is termed 
"open" cavity flow. Additional experimental studies of cavity-type 
separated flews are reported in refs. 5-10 for the laminar case and refs. 
11-15 for the turbulent case. 

Burgraff (ref. 16) approached the gap heating problem analytically 
for the laminar separation condition. His analysis postulated an invis- 
cid recirculating core flow within the cavity, and, therefore, is not 
applicable to fully viscous cavity flows. Nestler (ref. 17 and 18 ) 
approached the problem for both laminar and turbulent cases by analy- 
zing the shear layer development and reattachment, assuming a wall Jet 
heating decay from the reattachmenu point. 

The ,^ap heating problem for shuttle is one which will occur in 
deep gaps where the ratio of gap width-to-depth (w/d) is less than 0.08. 
(Gap width will be nominally - .082 cm.) Typically, the gap heating 
environment for shuttle will consist of an oncoming boundary layer with 
thickness far in excess of the local gap width (5*/v >> 1). All of the 
investigations previously mentioned, however, have dealt with cavities 
with width-to-depth ratios much greater than the shuttle value (usually 
> l), and boundary layer thicknesses generally much less than the gap 
widtn. Only the data of Weitirg (ref. 9) for laminar flow, begin to 
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approach the gap and boundary layer geometries of interest for the 
shuttle. 

In support of the space shuttle technology program, an experimen- 
tal effort has been focused on the shuttle-related gap heating problems. 
Evaluation of TPS tile material thermal performance (ref. 19) and defi- 
nition of the gap heating envii-onment are important elements in this 
effort. Johnson (ref. 20) studied the effects of gap width and boundary 
layer thickness on TPS gap heating for turbulent flow over a simulated 
tile array at Mach 8. This work revealed a potential heating problem 
area at the intersection of streamwise and transverse running gaps . 
Throckmorton (rel. 21) obtained data on a simulated tile array in a 
turbulent boundary layer to investigate the effects of gap width, flow 
angularity, and tile surface mismatch on tile gap heating. Foster 
et al (ref. 22) investigated the effect of similar parameters for single 
gaps and gap intersections in laminar flew. A compilation and analysis 
of these and other recent gap heating data may be found in ref. 23. 

Because of entry angles of attack and the complex curvatures of the 
orbiter external surface, boundary layer flows over much of the flight 
vehicle will be strongly influenced by gradients in surface pressure. 
Each of the recent studies of TPS gap heating phenomena, however, have 
been conducted on sharp flat plates or in tunnel wall boundary layers 
under conditions of zero pressure gradient. The present investigation 
was undertaken to assess possible effects of pressure gradient on the 
heat transfer to RSI tile array gaps under thick, turbulent boundary 


layer conditions. 
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IV. EXPERIMENTAL APPARATUS AND PROCEDURES 

Heat transfer and pressure measurement tests were made on a curved 
array of full-scale simulated RSI tiles submerged in a thick, turbulent, 
tunnel-sidewall boundary layer. Transverse pressure gradients of 
varying magnitude were induced over the model surface by rotating the 
curved array with respect to the flow direction. This enabled the study 
of RSI tile array gap heating as affected by pressure gradient. The 
tile array was tested both with gaps present, and with the gaps filled 
and smoothed to provide smooth-surface reference data. 

Heat transfer and pressure measurements were also made on a smooth 
flat plate mounted in the tunnel sidewall. These data, along with 
measured pitot pressure profiles, provided a definition of the charac- 
teristics of the 3-dimensional boundary layer flow in which tile array 
tests were conducted. 


Facility 

The experimental results presented herein were obtained in the 
Langley Research Center Continuous Flow Hypersonic Tunnel. This facil- 
ity, which has a 78.7*+ cm (31-inch) square test section, operates at a 
nominal 'f'reestream Mach number of 10.3 over a freestream unit Reynolds 
number range of 1.5 - 8.2 x 10^ per meter using air as the test gas. 

The tunnel may be operated in either a blowdown or continuous, closed- 
circuit mode. A schematic of the tunnel circuit is shown in figure 1. 
For continuous operation, the high pressure air supply and vacuum sphere 
are used to initiate the hypersonic flow. When the flow has been 
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established, the second minimum is closed down, and the compressors 
brought into the loop to maintain the hypersonic flow. Air heaters 
are of the tubular, electrical resistance type. The tunnel throat, 
expansion, and diffuser sections are all water cooled. A photograph of 
the facility is presented in figure 2. 

For these tests, the models were mounted on the model injection 
mechanism which is shown adjacent to the test section in figure 3- 
This device allows a model to be isolated from the hypersonic airstream 
for model cooling or geometric changes. The mechanism also provides 
rapid injection of & model into the hypersonic airstream. 

Models and Instrumentation 

Pressure gradient model .- The RSI tile pressure gradient model was 
a curved array of simulated full-scale RSI tiles fabricated of Inconel 
sheet with a nominal thickness of O.OU83 cm. Surface curvature was 
defined by a right circular cylinder of 102.28 cm radius, cut by a 
plane at a 5 degree angle to the axis of the cylinder. A model sche- 
matic is shown in figure h. Individual tile size was it. 92 cm square. 

The model was fabricated such that each transverse row of simula- 
ted tiles was a continuous sheet of material with the streamwise gaps 
formed by bending. The transverse gaps were fabricated separately and 
electron beam welded between tile row sections to form the complete 
tile array. Stress relief in the material following the welding pro- 
cess resulted in a "pinching down" of the transverse gaps from the 
desired width. Resulting gap width between tiles was nominally 0.30 cm 
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for streamvise running gaps and 0.20 cm for transverse gaps. Gap 
depth was 2.86 cm. A photograph of the model is shown in figure 5* 

Spacer plates were also fabricated to maintain model edges flush 
with the tunnel sidewall Three sets of spacer plates were fabricated 
to allow testing at model flow angularities of 0, ±22.5, and ±^5 degrees. 
The model and spacer plates mounted on the facility injection system 
for testing are shown in figure 6. 

The model was instrumented with 67 chromel-alumel thermocouples 
spot-welded to the model back-surface at tile surface and gap locations 
defined in figure 7- The model was also fitted with 12 static pressure 
orifices located as shown in figure 7. 

For smooth surface testing of the model, the gaps were elimina^d 
by filling them with twine to a point approximately one gap width below 
the tile exterior surface. The remainder of the gap was filled with 
plaster of paris which was allowed to dry and then smoothed by sanding 
to match the surrounding surface. This created a smooth model with 
correct curvature for measurement of surface reference data. 

Flat plate model .- The smooth flat plate model, used to measure 
undisturbed wall heat transfer coefficient, was fabricated from 321 
stainless steel sheet with a nominal thickness of 0.127 cm. The 50.8 cm 
square panel fit flush with the injection plate fabricated for testing 
of this model. The model and injection plate are shown mounted on the 
injection strut, ready for testing, in figure 8. The flat plate model 
was instrumented with 11 chromel-alumel thermocouples spot-welded to 
the model back-surface, located as shown in figure 9* 
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Boundary-layer probe .- A boundary-layer rake with 11 pitot probes 
was fabricated for use with the flat plate model. The rake could be 
placed downstream of the flat plate model at four spanwise locations. 
Individual pitot probes were located at stations normal to the wall 
as shown in the rake schematic of figure 10. Tube diameters varied with 
position relative to the wall. The rake is shown installed with the 
flat plate model for testing in figure 8. 

Test Procedures and Conditions 

The transient calorimeter technique was used to measure the heat 
transfer rates to the surfaces of the thin skin models. The tests 
were conducted with the models initially at room temperature, isolated 
from the hypersonic airstream within the injection chamber, at a pres- 
sure equal to the test section static pressure. With the hypersonic 
flow established in the test section, the model was rapidly injected to 
the test position, flush with the tunnel sidewall, and data were auto- 
matically recorded at a rate of 20 samples/second. After the model had 
been exposed tc the airstream for an interval of time sufficient to 
allow pressure transducer outputs to "settle out," the model was retrac- 
ted t m the stream into the injection box. 

Both the tile array and flat plate reference models were tested 

at nominal freestream total pressures of 2.1*1, 5.17, and 9.65 x 10^ 

2 

N/m ; corresponding to freestream unit Reynolds numbers of 1.6, 3-3, 
and 6.1 x 10^ per meter. The tile array model was tested at flow 
angles, a, of 0, ±22.5° , and ±1*5°* Testing at positive and negative 
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flow angles allowed measurement of both temperature and pressure data 
over the entire model surface while instrumenting only one side of the 
plane of symmetry for either temperature or pressure as shown in figure 
7. By "mirror imaging" the data obtained at negative flow angles, data 
for the entire model surface was obtained for flow angles of a - 0°, 
22.5°, and 45°. 


Measurement Techniques 

Temperature data.- Temperature data were obtained using the chromel- 
alumel thermocouples with a reference junction of 324.8 K. The inf- 
erence Junction temperature was thermostatically controlled to within 
±0.55° K. Freestream total temperature was measured at the centerline 
of the tunnel settling chamber immediately upstream of the nozzle. 

Pressure data .- Freestream total pressure was measured in the 

settling chamber by three strain-gage-type transducers. The transducers 

6 2 

had ranges of 0-3.45, 0-6.89, and 0-17.24 x 10 N/m . The transducer 
with the lowest pressure range which remained on scale was used to 
record the data. 

Static pressures were measured using baratrons mounted on the in- 
jection strut immediately behind the model. The baratron is a capacitance- 

type transducer, operated in conjunction with a signal conditioner to 

2 

allow measurement of pressure over seven ranges from 0-68.9 N/m through 
4 2 

0-6.89 x 10 N/m . The signal conditioner provides automatic ranging 
such that the measurements were obtained on the lowest possible range. 
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Data Reduction 

Freestream Flov Quantities .- Freest ream Mach number was determined 
from previously obLained tunnel flow calibrations. Freestream tempera- 
ture and pressure quantities were calculated using the one-dimensional 
perfect gas relations of ref. 2h with corrections for real gas imper- 
fections. The real gas correction factors were derived from the work 
of Erickson and Creekmore (ref. 25) on thermodynamic properties of 
equilibrium air. Fluid viscosity was calculated using the Sutherland 
relation for lew temperature, shown below in SI units: 


V = I.U 58 x 10 


-6 


^3/2 

X 

T + 98.33 


N-sec 

2 * 
m 


T - °K 


( 1 ) 


Boundary Layer Profiles .- Boundary layer velocity and density 
profiles were derived from flat plate wall boundary layer pitot pressure 
measurements. These calculations required assumptions of the distri- 
bution of static pressure and temperature through the boundary layer. 

As shown in figure 11, the measured flat plate static pressure was 
significantly higher than the computed freestream static pressure. 

Static pressure within the boundary layer was assumed to vary linearly 
between the measured wall and computed freestream values. 

Total temperature within the boundary layer was assumed to vary 
as the square of velocity, 



(2) 
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which is characteristic of turbulent nozzle-wall boundary layers 
(ref. 26). Unpublished total temperature measurements made by D. H. 

Crawford in the wall boundary layer of this facility show good agree- 
ment with the quadratic temperature-velocity relation. 

With the measured pitot pressure and assumed static pressure at 
each point, the Raleigh pitot equation was applied to calculate local 
Mach number. Local velocity and density were then computed using the 
assumed temperature distribution and the perfect gas equation of state. 

Heat Transfer Data .- The test procedure of rapid injection o** the 
isothermal model to the test position provided a step input in heat 
transfer to the thin-skin model. Heat transfer rates were determined 
by the transient-calorimeter technique of measuring the time-rate-of- 
change of the model skin temperature^ For data reduction purposes , the 
one-half second interval of temperature data immediately following mode* 
injection was disregarded to allow steady-state conditions to stabilize 
in the gap flow. This time is in excess of the required time as reported 
by Nicoll, ref. 27. A quadratic least squares curve was fit to the 
subsequent U-second interval of data for each thermocouple. Rates-of- 
change of temperature with time (3T /3t) were evaluated analytically 
from the curve fit expressions at the initial point of each curve fit. 

i * 

""A detailed assessment of the accuracy of this technique for these > 

tests is contained in Appendix A. 
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Heat transfer rates were then computed from the expression 

<1 = f ) m c p m X OV^t) (3) 

Heat transfer data are expressed in the form of the heat transfer 
coefficient (h) defined as 


h s f ' - T 00 

aw w 

Adiabatic wall temperature (T av ) was computed from the relation 



( 5 ) 


where recovery factor (r) was assumed equal to 0.89. 
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V. RESULTS AND DISCUSSION 
Coordinate Definition 

In the discussion of experimental results which follows, the data 
are referenced to coordinate sy- terns defined in figure 12. All flat 
plate and wall boundary layer data are referenced to a stream-oriented 
coordinate system (X, Y, Z) which is fixed within the flow. All curved 
panel data, both smooth and with gaps present, are referenced to a 
coordinate system (x, y, z) fixed within the rotating curved array. 

The origins of both coordinate systems are located on the tunnel side- 
wall at the center of rotation of the tile array model. 

Boundary Layer Surveys 

Boundary layer velocity profiles measured on the sidewall center- 
line are presented in figure 13 for the four freestream unit Reynolds 
numbers at which flat plate data were obtained. The profile shapes are 
characteristic of a fully developed turbulent boundary layer and the 
profiles tend to "fill-out" (boundary layer thins) with increasing unit 
Reynolds number. Profiles obtained at several transverse locations at 
a single unit Reynolds number are compared in figure ll*. A significant 
change in the boundary layer profile shape as a function of transverse 
position is observed. This transverse variation in boundary layer 
profile is an indication of the 3-dimensional character of the boundary 
layer flow in a nozzle of square cross-section. Although they do not 
present all of the measured pro^l^ data, figures 13 and lk illustrate 
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the trends observed in the data for the unit Reynolds numbers and 
transverse locations at which profiles were measured. 

In order to more readily define the effects of freestream unit 
Reynolds number and transverse location on boundary layer characteris- 
tics, the measured profile data were integrated to obtain values of the 
bulk quantities of boundary layer displacement (6 # ) and momentum (0) 
thicknesses. The quantities are defined by the relations 


00 



and are a measure of the mass flow and momentum flux deficits respec- 
tively, within the boundary layer. The variations of these quantities 
with freestream unit Reynolds number and transverse position are pre- 
sented in figure 15. Momentum thickness decreases, as expected, with 
increasing unit Reynolds number and shows little transverse variation. 
Boundary layer displacement thickness does not change significantly 
with unit Reynolds number. Displacement thickness does, however, 
decrease rapidly with transverse position, reflecting the changes in 
boundary layer profile shape noted in figure lU. 


Flat Plate Results 

Heat transfer data.- Measured flat plate heat transfer coefficients 


are shown in figure l6 for the full range of test conditions. At each 



1 


r 


1 


i 


i 


! 

15 

total pressure, the heat transfer increases in the transverse direction, 
primarily <^s a result of the thinning of the boundary layer, and is 
symmetric about the centerline. The streamwise distribution of heat 
transfer exhibits a heating decay characteristic of a fully-developed 
thickening boundary layer. 

Pressure data .- The flat plate pressure data (figure 17) show a 

transverse distribution similar to that of the heat transfer data. 

Streamwise, however, the pressure decreases to a minimum and increases 

in the downstream (+X) direction. This behavior is attributed to the 

fact that the test section area of this facility does not represent a 

complete expansion of the flow from the tunnel throat. The tunnel 

nozzle expands the flow to a point, approximately one meter upstream 

of the test section center, at which the flat, parallel walls of the 

test section begin. The intersection of the expansion section and test 

section walls then constitutes a compression corner for the wall boundary 

layer flow. The data presented herein were obtained in the compression 

region downstream of this comer. 

Flat plate data correlation .- As a basis for understanding the 

mechanisms controlling the flat plate wall boundary layer heat transfer, 

the heating data of figure 16 were correlated with parameters which could 

be expected to influence heat transfer. Figure 18 presents flat plate 

heating data expressed in the form of Stanton number, N , plotted as 

ot 

a function of Reynolds number based on wall conditions and boundary 

layer momentum thickness. Re ft . The heat transfer data points were 

v,o 
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inferred from the data of figure l6 for the transverse locations at 
which boundary layer profiles were obtained. 

At each transverse location, Stanton number and Reynolds number 
correlate the data over the range of freestream flow conditions. Tnis 
result is not unexpected, since for fully developed turbulent boundtuy 
layers over flat plates , heat transfer over a Reynolds number re 
can normally be correlated by Stanton and Reynolds numbers (ref. 

For this tunnel wall boundary layer flow, Stanton number was found 

to vary as the wall Reynolds number. Re fi , to the -.07 power at each 

w,u 

transverse location, as indicated in figure 18. Failure to correlate 
the data for the three transverse locations is not surprising, as the 
Stanton/Reynolds number correlation relates changes in streamwise 
variables and has no application to transverse flow phenomena. The 
upstream histories of the boundary layers affecting each transverse 
location are unique and no simple boundary layer parameter relates the 
transverse influence of boundary layer characteristics. The single 
parameter which varies significantly in the transverse plane is boundary 
layer displacement thickness; and for a given boundary layer flow, heat 
transfer is known to decrease with increasing boundary layer thickness. 
As previously noted, for this tunnel wall flow, displacement thickness 
does not vary significantly with unit Reynolds number. The effect of 
transverse boundary layer thickness variation on tunnel sidewall heat 
transfer is indicated in figure 19 where the data of figure 18 are 
presented as a function of displacement thickness. The flat plate heat 
transfer data correlate with pertinent flow variables according to the 
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relation 
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Re 


•07 = A «! 
w,e w 


+ B 


where W is defined as the tunnel test section 
half-width (39*37 cm). 


A = -0.0109 
B = O.COU66 


RSI Tile Array Results - Smooth Model 

Because of the complexity of the wall bo^mdary layer flow indi cated 
by the flat plate and boundary layer probe results , heat transfer ana 
pressure data were obtained on the RSI tile array model with no gaps 
present. The data were used to relate measured surface and gap heating 
to an undisturbed surface reference. 

Pressure data .- Smooth model pressure data and fourth-order least 
squares curve fits of the data are presented in figure 20 for all .est 
conditions. The data from the orifice located at y = ±5.08 cm were not 
considered for computation of the curve fit expressions as data from 
this transducer were consistently higher than that from the other trans- 
ducers. At an array rotation angle, a, of 0°, the pressures show the 
same transverse trends demonstrated by the flat plate data. Flow de- 
flection angle, $ (i.e. the angle between the freestream velocity vector 
and the plane of tangency Oj. the surface at a point), has negligible 
point to point variation between orifices when a = 0°. With ir.reases 
in rotation angle, a, to 22.5° and 45°, flow deflection angles on the 
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upstream portion of the model increase, those on the downstream portion 
decrease, and static pressures vary accordingly.^" 

Heat transfer data .- Measured smooth surface reference heat trans- 
fer data are presented in figure 21. For the zero rotation angle case, 
heat transfer increases, as expected, with transverse distance from the 
center? .ne. The rate of increase, however, is greater than that mea- 
sured for the flat plate. Hiis more rapid increase is hypothesized to 
be the result of an effective "thinning" of the boundary layer due to 
the protrusion of the model into the boundary layer flow. With a bound- 
ary layer which is much thicker than the protruding height of the tile 
array model, it is not thought that the. model will significantly affect 
the outer portions of the boundary layer as would be the case if the 
characteristic dimension of the model was of the same order as the 
boundary layer thickness. Rather, the boundary layer edge location 
remains essentially unchanged from the flat plate case and the boundary 
layer thickness is decreased by the protrusion of the curved model into 
the flow. 


^Note that the v-coordinate is fixed in the rotating array, while the 
Y-coordinate is fixed in the tunnel sidewall. While y-values are 
constant for each orifice, the orifice locations within the complex 
boundary layer flow vary with array rotation angle. Therefore, the 
pressure and heat transfer variations shown in the figures as functions 
of the y-coordinate, are a superposition of effects of changing flow 
deflection angle, and changes in Y-coordinate position (boundary layer 
conditions affecting the point). Derivation of the expression for 
flow deflection angle as a function of y-location and array rotation 
angle, a, is contained in Appr ix B. 



\ 




19 

Smootn surface heating data for the a * 22.5° and 45° cases show 
increased surface heating with increasing flow deflection angle 
(y-increasing) as expected. The data also show, however, increases in 
heating where the flow deflection angle is decreasing (y-decreaslng) . 
This anomaly is explained by the fact that these heating increases, in 
regions of decreasing flow deflection, are occurring at wall positions 
where boundary layer thickness is decreasing. The opposing effects 
of decreasing boundary layer thickness and flow deflection angle are 
dominated by the boundary layer thinning effect, and heating increases. 
In the following section, a method will be developed to predict smooth 
model surface heat transfer. This prediction method uses the measured 
flat plate heat transfer and pressure data as a base and perturbs these 
data to account for the flow deflection and boundary layer thinning 
effects which result f.om curved panel rotation. 

Prediction of smooth surface heating characteristics .- Consider 
the correlation of flat plate heating data presented in figure 19: 


W" * A V 


♦ B 


( 8 ) 


and assume that a correlation of this fora is valid for the smooth 
curved panel data. Therefore, for curved panel data: 
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where 6*^ is an effective boundary layer thickness as hypothesized 
in the previous section. By definition. 


St p u c 
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and Re 
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Substituting these expressions in Equation (9), 
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Assuming a perfect gas. 
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( 11 ) 


and assuming vail temperature fluctuations are negligible. 
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then. 
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Neglecting the veak dependence on 6: 



( 13 ) 


For the purpose of demonstrating the predicted distribution of 
heat transfer to all tiles in the model array, an estimate of the local 
pressure level must be made. Modified Newtonian Theory is applied to 
estimate the wall static pressure variation as affected by transverse 
location and flow deflection angle. The Modified Newtonian expression 
for pressure, referenced to the flat plate value (0 = 0) is: 



Now ^timate 


(lM 


6* = S* - kz _ 

eff surf 


( 15 ) 


where z . is the local protrusion of tae smooth model surface into 
surf 

the boundary layer flow and k is a constant. The value of ohe con- 
stant, k, which provided a "best fit" to the experimental data was 
found to be 0.5. 

By applying the known flat plate pressure, heat transfer, and 
boundary layer information to expressions ( 15 ), (lM, and ( 13 ), the 
distribution of heating to the surface of the smooth, curved panel may 
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be predicted for all flow conditions and panel rotation angles. Pre- 
dicted heat transfer distributions for the smooth model at a freestream 
Reynolds number of 3*3 x 10^ per meter and array rotation angles of 0° , 
22.5°, and 45° are shown in figure 22. Comparison of the theory of 
figure 22 with the corresponding data of figure 21 demonstrates that 
the prediction of transverse heating trends by Equuti (13) is good. 

The failure of the theory to more accurately predict heating level is 
attributed tc the inaccuracy of the modified Newtonian pressure predictions. 

A comparison is made in figure 23 of the measured and predicted 
smooth panel pressures. The failure of the modified Newtonian method 
to more accurately predict pressure level was not unexpected as this 
method is applicable to hypersonic flows of uniform total pressure. In 
this case, the method was applied to a boundary layer, a flow of non- 
uniform total pressure, without consideration of the total pressure 
variation within the boundary layer. The intention here, however, was 
to reproduce trends in pressure level, not necessarily magnitude; and 
the modified Newtonian approach did this satisfactorily. 

If the actual local pressure level were known, it then appears 

that the local heat transfer could be predicted reasonably well using 

this approach. Figure 24 presents measured heat transfer data for the 

transverse ray along which pressure data were also measured for the 

3.3 x 10** /m Reynolds number case. The predicted heating from Equation 

(13), using the measured pressure data, shows excellent agreement with 

the measured heat transfer. Similar results were observed in the data 

, ^ 6 

obtained at Reynolds numbers of 1.6 and 6.1 x 10 per meter. 
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RSI Tile Array Results - Gaps Present 

Pressure data .- Measured pressures for the simulated RSI tile 
array are presented in figure 25 for all test conditions. Tile surface 
data are shown as open symbols, gap floor data as solid symbols. Com- 
parison of these data with that of figure 20 indicates that the surface 
pressure distributions with gaps present are essentially identical to 
those obtained with a smooth model. Pressure level within the gaps 
closely follows that of the tile exterior surface. 

Gap heat transfer data .- Measured distributions of heat transfer^" 
along the instrumented transverse gap are shown in figure 26 for each 
array rotation angle at the 3.3 x 10^/m Reynolds number flow condition. 
For the zero rotation case (figure 26(a)), at a transverse position off 
the centerline, the surface flow completely bridges the gap with no 
intense reattachment heating occurring at the gap corner, nor any sig- 
nificant change in surface heating downstream of the gap. Note also 
that off the centerline at a depth into the gap of lesB than four gap 
widths (^-symbol), the local heating level is less than 2 percent of the 
undisturbed surface value. The solid symbols of figure 26, and those 
figures to follow, indicate heat transfer data which are of questionnable 
accuracy. The model skin temperature increases from which these data 
were derived were of the same order as the precision of the measurement 
apparatus. The heat transfer rates which resulted from the temperature 

^All heat transfer data discussed in this section are non-dimensionalized 
by the measured heat transfer to the smooth surface model. Surface data 
are normalized by the smooth model measurement of the same thermocouple; 
gap wall data are normalized by the smooth model surface measurement at 
the thermocouple location nearest the gap. 
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data curve fits are then affected as much by recording system data 
scatter as by steady state con/ective heating. These data are included 
for completeness. 

For the zero rotation case, figure 26(a), on the array centerline 
(y s 0), the gap geometry consists of a streamwise gap intersecting the 
transverse gap, creating a stagnation region on the transverse gap wall 
at the intersection point. Heat transfer within the gap at the inter- 
section point is significantly higher than that for a simple transverse 
gap as exists away from the centerline. This heating increase at the 
gap intersection is attributed to impingement of the "channel" flow in 
the streamwise gap upon the forward-facing wall of the transverse gap. 
Excessive heating at such gap intersections poses a serious problem for 
the shuttle TPS designer. Dunavant and Throckmorton (ref. 29), using 
gap intersection data from several facilities , have shown that such 
data may be correlated as a function of boundary layer displacement 
thickness, streamwise gap running length, gap width, and gap depth 
location. 

Rotation of the tile array, figures 26(b-c), produces slight 
increases in the surface reattachment heating, as effective gap width 
increases with rotation, and negligibie changes in the heating levels 
within the gap. The only significant change in the gap heating dis- 
tribution is a transverse shift of the heating peak within the gap which 
results from the gap intersection phenomena , as array rotation angle 
increases. Results obtained at Reynolds number of 1.6 and 6.1 x 10^/m 


are similar to those shown. 
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Figure 27 presents heat transfer data for the streamwise gap at 
a Reynolds number of 3.3 x 10^/m as a function of array rotation angle. 
The significant variations in gap heating observed in these plots result 
from the transition of the wall from one which is upstream-facing to 
one which is downstream- facing as array rotation angle changes. When 
the wall is downstream-facing (a > 0), the flow appears to separate up- 
stream of the corner thermocouple, resulting in heating lover than the 
undisturbed surface value. When the wall is upstream-facing (a < 0), 
the flow reattaches in the vicinity of the corner thermocouple resulting 
in heating equivalent to or greater than the undisturbed surface value. 
This result is in opposition to that observed in figure 26(c), (y = 3.81 
cm), where the flow appeared to completely bridge the gap, attaching 
downstream of the corner. These contrasting results are attributed 
to differences in the gap width and corner radii between the streamwise 
and transverse gaps. As a result of the model fabrication process dis- 
cussed previously, streamwise gap width was approximately 1.5 times as 
large as transverse gap width and streamwise gap edge radius was sig- 
nificantly larger than that of the transverse gap. The increased gap 
width and large exposed surface area at the streamwise gap edge allow 
for diffusion of the shear layer into the gap and flow reattachment at 
the tile corner as opposed to the "bridging" phenomena observed for the 
transverse gap. 

The variations of gap heat transfer with rotation angle for the 
simple transverse gap, and the gap intersection, are presented in 
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figures 28 and 29, respectively, for the 3.3 x 10^ /m Reynolds number 
case. Array rotation angle has negligible effect upon the heating to 
the simple transverse gap (figure 28). However, for the gap intersec- 
tion location (figure 29), array rotation tends to increase the local 
heating near the tile surface while resulting in decreased heating at 
depth within the gap. Again, data at other unit Reynolds numbers 
demonstrate similar trends. 

Pressure Gradient Effects on Gap Heat Transfer 

The discussion of the previous section has centered upon the 
effects of array rotation on the heat transfer within the gap. 

Rotation of the curved array not only produced crossflow over the 
array, but also the desired pressure gradient within the gap, providing 
for determination of the effect of pressure gradient upon heat 
transfer within the gap. The magnitude of this pressure gradient was 
determined by differentiating the curve fits of the gap pressure data 
of figure 25. 

In figure 30, measured heat transfer in the transverse gap at 
y = ±3.8l and ±11. U3 cm is presented as a function of a non-dimensional 
pressure gradient parameter. This parameter, (3P/3y) (L/P), is the 
local pressure gradient divided by the local pressure times some charac- 
teristic length. Using the full-scale tile dimension as the character- 
istic length (L ■ 15. 24 cm), this parameter is physicalxy a AP/P 
parameter for one tile length along :he gap. No systematic effect of 
pressure gradient on gap heat transfer is evident. 
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Shuttle Pressure Gradient Simulation 
In order to evaluate the full-scale shuttle pressure gradient 
simulation obtained with the curved array model, full-scale vehicle 
spanvise pressure distributions were obtained at two longitudinal sta- 
tions. These pressure distributions were calculated for a full-scale 
vehicle at Mach 10.0 and 30° angle of attack in a real gas. The surface 
pressure distributions were numerically differentiated to obtain the 
full-scale values of the pressure gradient parameter, (3P/3S) (L/P). 
Values of this parameter for the full-scale vehicle, are presented in 

figure 31 as a function of surface dimension, S/Y , for the two longi- 

max 

tudinal stations illustrated. Hie shaded area superimposed upon these 
data indicates the range of pressure gradient parameter values obtained 
in the wind-tunnel tests. Simulation of full-scale vehicle pressure 
gradient levels is excellent with the exception of the wing leading 
edge regions. Ving leading edge TPS, however, will be a solid material 
with no gaps present. Therefore, the effect of pressure gradient on 
gap heating is not of concern on that portion of the vehicle surface. 
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VI. CONCLUDING REMARKS 

An experimental investigation has been undertaken to assess the 
effect of pressure gradient on the heat transfer to reusable surface 
insulation (RSI) tile array gaps submerged in a thick, turbulent bound- 
ary layer. The experimental program consisted of heat transfer and 
pressure measurements on a curved array of full-scale simulated RSI 
tiles in the tunnel wall boundary layer of the Langley Research Center 
Continuous Flow Hypersonic Tunnel over a range of freestream Reynolds 
numbers and flow angularities. The tile array model was tested with gaps 
present, and with gaps eliminated to obtain smooth surface reference 
data. In order to gain a thorough understanding of the boundary layer 
flow in which these tests were conducted, measurements of tunnel wall 
boundary layer pitot profiles and flat plate wall pressure and heat 
transfer were also made. 

The measured flat plate wall heat transfer data were correlated 
with measured boundary layer parameters and a method was developed for 
prediction of smooth, curved array surface heat transfer. An assess- 
ment was made of the full-scale vehicle pressure gradient simulation. 

The results of this investigation indicate the following: 

(l) The wall boundary layer in the square test section of the 
Langley Research Center Continuous Flow Hypersonic Tunnel is highly 
three-dimensioned with significant transverse thickness variations. 

{2) Heat transfer to the tunnel wall correlates with momentum 
thickness Reynolds number and boundary layer displacement thickness. 
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(3) Heating to the curved tile array in this flow was more sensi- 
tive to transverse position than was the flat plate data. This sensi- 
tivity resulted from thinning of the boundary layer by the protrus* 'in 
of the curved array surface into the flow. 

(b) A method has been derived for prediction of smooth, curved 
array heat transfer in this tunnel wall boundary layer flow. The 
method uses predicted or measured surface pressure to perturb the mea- 
sured flat plate heating data, and accounts for "effective" boundary 
layer thinning due to the protrusion of the model into the boundary 
layer. 

(5) The level and distribution of pressure on the floor of the RSI 
tile array gaps follow closely those of the external tile surface. 

(6) Heat transfer to the gap wall is significantly higher at the 
intersection of a streamwise and transverse gap than for a single 
streamwise or transverse gap. 

(7) Simulation of full-scale vehicle pressure gradient was good. 

(8) For a thick turbulent boundary layer, there is no systematic 
effect of pressure gradient on tile array gap heat transfer. 
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Figure 1.- Schematic diagram of t facility. 
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Figure 8,- Flat plate model mounted on injection strut. 
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Figure 9*- Flat plate model instrumentation locations. 
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Figure 10.- Boundary layer pitot probe schematic. 
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Figure 11.- Comparison of centerline wall and free? cream statin pressures. 






Figure 12.- Definition of coordinate systems. 











L6 -8 0 8 16 2< 

Y, cm 


Figure 16.- Heat transfer coefficient distribution on the flat plate 























R< 







"St^w.O 



6 # 

W 


Figure 19.- Correlation of flat plate heat transfer with displacement 
thickness . 
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Figure 29.- Variation with array rotation angle of heat transfer 
to the transverse gap at its intersection with the 
streamwise gap. Re_/m = 3.3 x 10^. 




Figure 30.- Effect of pres 
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VIII. APPENDIX A 


Heat Transfer Rate Measurement Technique 


Aerodynamic heat transfer rate measurements were made using the 
thin-skin or transient-calorimeter technique. This technique uses the 
measured temperature response at the backface of a "thin" metallic 
model and the assumption that the model acts as a heat sink only with 
no front-to-backface temperature gradients (zero thickness) to 
compute model frontface heat transfer rates. Heat transfer rates 
computed by this method may be in error, however, because of the 
finite thickness of the model skin. The magnitude of this error may 
be assessed by comparing the solution for transient 1-dimensional heat 
conduction through a finite slab to the result for the "thin-skin" 
assumption. 


Finite Thickness Slab 
Consider the 1-dimensional transient 
conduction of heat through a finite 
solid of thickness X, density p, 
specific heat c, and thermal 
conductivity k: 




where k * 

pc 
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JL 
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Initial condition: 

T (x,0) = T. 

Boundary conditions: 

(Backface adiabatic) 
r Constant 


k 9T(0.t) _ Q 
* 3x 0 


3T(X.t) = 


3x 


( A— 3 ) 


(A-U) 

(A-5) 


The exact solution to this equation may be written (Ref. 30): 


T(x,t) - Ti - £ (terfc ISLlji 

n=0 


(2n + l)A - x 

ft* 


. ierfc (2n*lYh*_x | 

2 fit 

(A-6 ) 


where ierfc q = J* erfc S d| 


Solving for the backface temperature (x * 0): 

{ierfc 


n=0 


(A-7) 


Define 


_ Kt 
A2 


F = — 
* O — 


then 


T(0,t) - Ti = 



(2n + 1) 


(A-8) 
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Thin Skin Approximation 

The "thin skin" approximation assumes a: - infinitely thin wall 
and therefore a simple energy balance between input heating rate and 
heat storage: 

pel = q = constant (A-9) 

Initial condition: 

T(0) = Ti (A-10) 

The solution is then: 

T(t) - Ti = q ~ (A-ll) 

or 

T(t) - ^ ^ F 0 (A-x£; 

Comparison of equations (A— 8 ) and (A-12) yields a direct relation 
for "thin-skin" versus finite-skin temperature response at the model 
skin backface. The ratio of actual finite-skin heating rate (<l finite ) 
to "thin-skin" heating rate ( Q.-thin ) is a Action only of the 

parameter F 0 = and varies as shown in the following graph. 

A* 
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IX. APPENDIX B 
Surface Flow Deflection Angle 

z 



Consider the unit surface normal vector, n: 


/N 


n 



A 





and the unit freestream velocity vector, u: 


(B-l) 


U a (u x , u y , (B-2) 

The scalar product of tfc »sf. vectors defines the angle between them, 

fi*u «(fiu + hu + fifl)a cos (90 + $) 
x x y y 2 2 T 


(B-3) 
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where (90 ♦ <J>) is the angle between the freestream velocity and 
surface normal vectors, and 4> is the angle between the surface tangent 
plane and the freestream velocity vector. 

A A 

Resolving unit vectors u ana n into their component parts: 


u 

X 


-cos 

a 

u 

y 


-sin 

a 

z 

= 

0 



(B-U) 


and 





2 



\ 



1 


r 
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Substituting the components in (A-3) : 


i 


A A 

n*u = -cos 


8 cos a sin 5° - sin 8 sin a = ccs (90 + 4>) (B-6) 


But, 


'surf , __ „ _ y surf ^ 

sm 3 — 7T 1 " — and cos p — |1 — 


R 


surf 


. L . y J^L-\ 

V P S urf / 


( B— 7 ) 


Therefore, 


cos (90 + 0) = 


i'-’fty 

surf 


(B-8) 


or 


( ■ HT 


y V-2 

sin <J> = |1 SUI "? | cos o sin 5° + 


/ y surf\ „ 

VCT““ 


(E-9) 



